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Abstract. On compact manifolds which are not simply connected, we prove 
the existence of "fake" solutions to the optimal transportion problem. These 
maps preserve volume and arise as the exponential of a closed 1-form, hence 
appear geometrically like optimal transport maps. The set of such solutions 
forms a manifold with dimension given by the first Betti number of the man- 
ifold. In the process, we prove a Hodge-Helmholtz decomposition for vector 
fields. The ideas are motivated by the analogies between special Lagrangian 
submanifolds and solutions to optimal transport problems. 

1. Introduction 

In |KMW| , graphs of solutions to the optimal transportation were shown to solve 
a volume maximization problem, using a calibration argument. With the appropri- 
ate metric on the product space M x M, maximality is equivalent to the vanishing 
of certain differential forms along the graph of the optimal map. In this note, 
we discuss the converse and see that, at least in the smooth case, topology allows 
for maximizers of the volume problem which do not arise as solutions to the opti- 
mal transportation problem, despite locally having the same geometric properties. 
These maximizers are special Lagrangian in the sense of Hitchin IT and Mealy 
[Mej , a pseudo Riemannian analogue of the special Lagrangian geometry of Harvey 
and Lawson [HL| . 

In the case when the cost is given by Riemannian distance squared, Delanoe 
[D] introduced the notion of "Lie solutions of Riemannian transport equations" 
(see also [L].) The graphs of Delanoe's solutions are maximizers of the volume 
maximization problem discussed in }KMWj . 

We recaU the McLean Theorem |McH Theorem 3.6] |Ma| Theorem 3.21] 

Theorem 1.1 (McClean). Suppose L is an smooth embedded special Lagrangian 
submanifold of a Calabi- Yau manifold. The moduli space M of special Lagrangian 
submanifolds near L is a manifold of dimension bi{L). The tangent space to A4 is 
identified with the harmonic 1-forms on L, which has a naturally induced L^ metric. 

The seminal paper of Harvey and Lawson shows that special Lagrangian sub- 
manifolds are minimal submanifolds in Calabi- Yau manifolds, and that a manifold 
is special Lagrangian if and only if the Kahler form and a certain n-form vanishes 
along the submanifold. Hitchin [H] analyzed the metric on the moduli space in 
McLean's theorem and showed that it arises via a Lagrangian embedding into a 
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pseudo-Riemannian space. Hitchin's notion of special Lagrangian, that the Kahler 
form and certain combinations of n-forms vanish, describes Lie solutions of the 
mass transport problems. 

The optimal transportation problem is the following. Given probability volume 
forms p and p on manifolds M and M, and a cost function c : M x M M, find a 
map T : M ^ M which minimizes a cost integral, 



c{x, T{x))dp 

IM 

among all maps T which preserve the volume, i.e. 

(1) T,p = p. 

The work of Brenier [B] and McCann |McC] shows that given standard conditions 
on the cost function, the unique solution will be the map satisfying ([T|) and arising 
as the cost exponential of the gradient of a potential function u : 

(2) T{x) — c-exp,-^du{x). 
A local version of ([1]), namely, 

(3) T*p = p 
is equivalent to the form 

pdx — pdx 

vanishing on the graph {x, T{x)) C M x M, and ^ implies the vanishing of a 
Kahler form defined in .KM]- The more general problem we attack here is to find 
maps which locally solve the optimal transport problem, that is, satisfy ([3]) and 
arise as the exponential of a closed form 

(4) T(a;) = c-exp,(r;(a;)), dyy = 0. 



The result is the following. 

Theorem 1.2. Suppose M, M are compact manifolds with nowhere vanishing smooth 
densities p and p. Let c he a continuous cost function on M x M, which away 
from a cut locus C, is smooth and satisfies local and global twist assumptions (A2) 
and (Al). If a diffeomorphism T is a smooth Lie solution of the mass transport 
equation (i.e. satisfying and ^) avoiding the cut locus, then there is a mod- 
uli space of smooth Lie solutions near T which is a smooth manifold of dimension 
bi{M). 

Remark: There are cost functions such that even for some smooth densities, the 
optimal solutions will not be smooth, as was demonstrated by Loeper [Lo], when 
the (A3) assumption on the cost is not satisfied. Also, in the absence of the 
(A3) assumption, local c-convexity does not imply global, so even on simply con- 
nected domains. Lie solutions are not necessarily solutions to the optimal transport 
problem. 

The core of our result is a Hodge-Helmholtz type decomposition holding at Lie 
solutions, which decomposes deformations of maps into those which preserve prop- 
erty ([3]), those which preserve property ([4]) and the harmonic deformations, which 
preserve both. The metric used is the one induced by the linearized operator to the 
optimal transport equation, multiplied by an explicit conformal factor. This metric 
and concomitant Laplacian are defined in section 2, and their Hodge-Helmholtz 
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properties are shown in section 3. In section 4 we give a possible politico-economic 
application of solutions. 

Acknowledgment 1.3. The author is very grateful for advice from Phillipe De- 
lanoe. 

2. Preliminaries 

For n-dimensional manifolds M and M, let c : N C M x M Hhe sl continuous 
cost function which is smooth almost everywhere, except on a set C=M x M — N 
which we call the "cut locus." (The reason for this terminology is clear if we use 
the distance squared function on a manifold as the cost function.) Let DDc be 
the n X n matrix given by (in local coordinates) 

(DDc) (x, x) — TT—^ — -c(x, x). 

\ ^ I dx'-dxi ^ ' 

On N we will require that 

(A2) det {DDc) ^ 

which is a local version of the standard twist condition: For each x and x respec- 
tively, 

(Al) x^Dc{x,x) 

is invertible, with an inverse depending continuously on x. 

To be clear with our conventions, we recall the following Kantorovich problem: 

If 



J{u,v) = / {—u)dp+ / vdp, 
Jm Jm 

the problem is to maximize J over all —u{x) + v{x) < c{x, x). One also considers 
a dual problem: If 



/(tt) = / c{x,y)dTT, 
Jmxm 

find the minimum of / over all measures tt on the product space M x M which have 
marginals p and p. It is well known (cf fV) that 

sup J{u,v) — inf /(tt). 

~il{x)+v{x)<c{x,x) 'ren(p,p) 

With this setup in mind we can derive the optimal map T from u as follows: 
Suppose {xq,xq) is a point where the equality ~u{xq) + v{xq) — c{xq,xq) occurs. 
The function 

Zxo{x) = c{x,X(i) +u{x) 

must have a minimum at xq. Then define the cost exponential T{xo,du) — xq. If 
differentiable, from the fact that z^o^x) is at a minimum we have 

(5) Ui + Ci{x,T{x,du)) = 

where Ci refers to differentiation in the first variable. (One can check that T{x,u) = 
— expj, Vit when c(x,x) = d^{x,x)/2). This only depends locally on the function 
u, and clearly requires that Du stay inside the range of Dc{x, ■). The elliptic 
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optimal transportation equation can be derived by taking another derivative and 
then a determinant: 

(6) + Cij + CigTJ = 

(7) det(u,j + c,j) = det(-c,,-T/) = dct(-c,g) 

using the fact that T locaUy pushes p forward to p, so satisfies 

(8) p{T{x))AeiDT = p{x). 

Here and in the sequel we use the following conventions: Indices i, j, k etc will be 
coordinates on M while p, s, v, etc will be coordinates on M.The variable t will be 
reserved for variations. We assume on any local product chart that c^j is negative 
definite, and let bis — "C^j^ with fe^* its inverse. The tangent space is denoted 
TM. We use Wijlx) = Uij{x) + Cij{x,T{x)). To be clear, we will use -^Wij to 
denote derivatives of w, otherwise a subscript will denote differentiation. We note 
that most of the quantities we deal with are coordinate invariant. In particular 
Wij actually is an honest tensor. For future use, we make note of the following 
restatement of the above identity © 

(9) Wij = bigT. 



2.1. Linearizing the elliptic equation. In order to linearize ([7]) take a variation, 
(10) u{x)+tv{x). 

First, we insert (fTO|) into ([5]) and differentiate to obtain 

v,^b^-s{x,T{x,du))DtT'' 

hence 

DtT'' = T/ = v,¥\ 

Now, taking a logarithm of ([7]) 

F{x,Du,D'^u) = 

lndet(uij + Cij{x,T{x)) — In dot (6^5(2;, 7" (a;)) ~ lnp(a;) + lnp(T(x)) 

which is linearized 
d 

di' 

= w'^v - w''b,^g¥^Vk - bHg.pbP'^Vk + (Inp), b^'vk 

A version of this linearized operator was introduced by Trudinger and Wang |TW[ 
2.18]. 

2.2. The KM and modified KM metrics and a related Laplace-Beltrami. 

In (KM] Kim and McCann considered the following pseudo metric on the product 
space M X M : 

, / b,s 

and symplectic form 

ui = bisdx^ A c?a;^. 



(11) Lv = -F{u + tv) = w'^ (v ,j + c,,-sT!) - bHs^pTf + (In p), Tf 
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In this metric, the graph {x,T{x)) of the cost exponential of any locally c-convex 
potential u is space-like and Lagrangian. The induced metric, in terms of coordi- 
nates on M, is given by 

{x,T{x))*h = Uij +Cij = Wij. 
For given mass densities p and p, consider the following metric in |KMW) 

1 f PP V^" ( fo,. 



2 V det 6,., ) \ b% 
and the calibrating form 

The graph of any solution to the optimal transportation problem will be a calibrated 
maximal Lagrangian surface with respect to this metric, but the converse is not true 
in general. Calibrated maximal Lagrangian surfaces are what we are studying in 
this paper. 

To begin, we define another metric. We will show that with respect to this 
metric, the tangent space of deformations of calibrated submanifolds coincide with 
harmonic 1-forms, which is the same situation that occurs in McLean's theorem. 
The metric we use is yet another conformal factor of the metric used by Kim and 
McCann, differing by a power of the conformal factor from the metric associated to 
the calibration in |KMW) . For a given frame on M, and a solution of the problem 
T define 



(12) gij{x) ^ Wij{x) 



f p{x)-p{nx)) y/^--'^ 

\det hrs{x,T{x)) 



This metric appears naturalyl suited for our problem, as we will use the Hodge 
decomposition with respect to this metric in the sequel. 

Before proving the next claim, note that at any point x G M, one has by the twist 
condition that for any covector 77 e '^x^^ there is at most one point x = T{x,r]) 
in M so that 77 = -~D(x,T{x,r])) This defines T whenever 77 is in an appropriate 
domain. 

Our first claim is the following. 

Proposition 2.1. Let n > 3. Let T = T{x,ri) be a cost exponential of a 1-form 
locally given by du. Define 

6 — hi det Wij — In p{x) + In p(x, T(x)) — In det big, 

and 

p{x)-p{T{x)) 



^det bis{x,T{x)) 
then with the above metric U3\) on M, 



Lz = A(^AgZ + i(V0,Vz)g 

where AgZ is Laplace- Beltrami operator with respect to g and L is a the linearized 
operator defined by ill]) . 
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Proof. Defining 
with 



A = 



p{x)-p{T{x)) 



^det bis{x,T{x)) 
we compute in local coordinates: 

- T;wr^=(^"^'"'ydit^u;^^-.,). 
A"/-^-y/det Wij 

= \w'^z,j + \{-w'''w^^ + ^^-^w'^' (InA) - + iu;'='^Wfciw'J)z,. 

A A ax^ 2 ■'2 ax-? 

Locally, because r/ = du, we can write Wab = + Cab, thus 

AAgZ = w'^ z^j^~w'''w^^-^wki~-^wij)-w'''w'' -^wij + ^^-^w (In A)^- + iw''' 

= w'^z,, + S^w'Hkisb'' - bi,,¥'w'^ - ^w'^ (Indetu;)^. + ^^w'^ (lnA)^-| z, 
(13) 

w'^^bkisb^^ - bijsb^'w'-^ 



W Zi 



= Lz--X{V9,Vz)g 



1 

2' 

where we have used ^ repeatedly. In (IT51) we combined the following two identities 

In A = - — {lnp(a::) + lnp(a;, T{x)) - Indet 
(n-2) 

In det w = 9 + hi p{x) — In p{x, T(x)) + In det bis 
before taking the derivative. □ 

3. Deformations 

Let T : M — >■ M be a Lie solution of the transport equation. As was used in the 
previous calculation, this means that T is locally a solution of the optimal transport 
equation, or equivalently, T is a cost exponential of a closed 1 form, satisfying 

detDT- 



p{T{x)) 

Throughout this section we will use the metric g on M defined by ((T2|) . 
Given a 1-form on M we can define a (vertical) deformation vector field 

where ds is the coordinate tangent frame for TM (to be precise, F is a section of 
the puUback bundle T*TM). For a vector field in this bundle, define a deformation 
of the map T via 

Tv(x)=exp*^)(T/(x)). 
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(We have fixed arbitrarily a metric on M in order to define exp*^, which is no 
more than a convenience.) 

Along N = MxM — Cwe define the 1-form a = Ci{x, x)dx'^ (which differs from 
the total differential of c by Cs(x, x)dx'^ ). Then we get the following exact Kahler 
form 

oj = da = —Cis{x, x)dx^ A dx^ 

on N. 

At a solution T we define the map 

$ : Ai(Af) A^{M) © A°(M) 

via 

^'(77) = {{Id X Tv)*Lj, * {Tv{x)*pdx - pdx)) 

Note that the level set $~^(0,0) consists of forms whose corresponding maps 
Tv are Lie solutions. Equivalently, the map {Id x Ty) is a calibrated Lagrangian 
submanifold of M x M with respect to the metric in [KMWj . 

Lemma 3.1. The image of ^ lies in exact 2-for'ms and coexact 0-forms. 
Proof. The first factor is a pullback of an exact form. For the second factor, 

J Tv{x)*pdx — pdx — 0, 

which follows easily from the fact that the diffeomorphism T is simply a change of 
integration variables. We are assuming the total mass of both densities is 1. □ 

Lemma 3.2. AtQ<E C''+^'°'{A^{M)), the differential satisfies 

D^{ri) = {dri,d*r]). 
In particular, is a topological linear isomorphism 

d*C''+2'"(A2(M)) ©dC'=+^'"(A°(Af)) — ^dC'^+i'"(Ai(M)) ©d*C'=+i'"(Ai(M)). 

Proof. To compute the derivative, for the first factor of the map $ we use the Lie 
derivative on the manifold M x M (cf [Mai section 2.2.2.]) 

j^{Id X TvYujlt^o = C(^^c^^^^d,){Id X ry)*w 

= {Id X Tv)* [-c^"r?™95jdw - d{c'"'r,^ds^uj)] 
= {Id X TvT [d{-c'"'f]^dsju:)] 
= iIdxTvr[d{r]^dx"')] 
— drj. 

For the second factor, to begin we note that because the volume element is given 

by 

/ / - \ 11/71 — 2 / / _\n/n — 2 

we have 

*{T*pdx - pdx) = ^ (det DTvp{T) ~ p) . 
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Differentiating 
1 



^ deiDTvp{T)-- 
pX A 

(14) = dctDTvp(T)^{(lndetDTy + (\np)' - (Inp)' - InA'j + ^ (In A)' . 

p\ A 



Noting that at 



det DT = 



piT) 



the expression p4)) becomes 

\ {(lndetOT)' + (hip)' - (Inp)'} . 
A 

Now in particular, differentiating with respect to t we have 

^ * {T^pdx - pdx)\t^o = \ \{DT-%Tl, + (lnp),T, 



[DT-X{b''^r,k)^ + {\np),b''vk 



= J [b.-sw^^b^'r^kj - h,w'^lfH'''-{b^p,+b^pr¥'wij)m + (lnp),6^"' 

(15) = ^ [w^'m^j - w'^bPH,p,7jk - bP%prb"ijk + (Inp), b^'rjk] 

(16) = d*77 

using ([5]) repeatedly. The last equality follows from the proof of Proposition 12.11 

Xd*r, = w''ij,,j + {-kjsb'W^ + (lnp),6"- - b'%spb^'} n,. 

The latter conclusion of the Lemma follows from standard algebra using the Hodge 
decomposition 



C''+^'"{A\M)) = ©rfC'=+2'"(A°(M)) ©d*C'^+2'"(A2(M)). 



□ 



Proof of Main Theorem when n > 3 

smooth map $ has surjective differential 



Lemma 



shows that at 0, the 



. C''=+i'"(Ai(M)) d*C''+^'"{A^{M))(SdC''+^-°'{A\M)) 

onto the product of exact and coexact forms. The kernel at is the harmonic 
forms, which splits by the Hodge decomposition. It follows by the Implicit Function 
Theorem (see [Ma, Thm 2.11]) that for each harmonic form 77 close enough there 
is a unique form xiv) lyi^ig in the orthogonal complement of the harmonic forms 
so that $(77 + xiv)) = 0- Thus a neighborhood of in the harmonic forms on M 
parametrizes the moduli space near T. 
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3.1. n = 2.. When n = 2 we are missing the power of the Hodge theorem, which 
asserts that the kernel of the operator has the same dimension as the first coho- 
mology class. We must somehow show that the set of solutions of { ([T5)l = 0} has 
this same dimension. Forming the elliptic operator 

L^d*d+S*S 

where 5 : A-'^ — A° is the operator defined by p5|) . the Implicit Function Theorem 
arguments as used above provide that the dimension of the moduli space near a 
solution is the dimension of the kernel of L. 

For a given , consider the transportation problem of finding 

(17) F {x, 9) = (T (x, d) , e {x, 9)) : x ^ x 
minimizing the cost 

c[(x, 9) , {x, 9)] = c{x, x) + distil (9, 9) 

among all maps (F, Q) pushing p A d9 forward to p A d9. 
For any Lie solution T : ]VP ]VP it is clear that 

(18) Fix,9)^iT{x),9 + 9o) 

is a Lie solution to the problem ([TTl) . By Theorem 11.21 for n = 3 there is a space 
of deformations of solutions, and it has dimension equal to 1 + bi{M). Suffice then 
to show that these deformations decompose into deformations in each factor. 

Let T be a Lie solution on M and let (IT51) be a Lie solution on M x S^. Let tj be 
a 1-form which defines a tangent vector to the space of deformations on x 5'-'^ , 
at the solution F, and write 

77 = r/iUJ^ + ri2Uj'^ + r]0d9 
for some cotangent frame w^jW^ for M^. Differentiating the equation 

d*r] = 

in the 9 direction, we have that 

w''v^.Je + {-bb,sb''w'' + (Inp), b'' - b'%kpb^'} r,,,e = 0, 

using the fact that the warped product metric does not depend on 9. Also, because 
77 is a closed form, locally we have rjg^ij = rjijg and r]i,g = ijg^i. Thus (the honest 
function) z — rjg locally satisfies an elliptic equation of the form 

w^^ Zij + A'zi = 

so enjoys a maximum principle. We conclude that on the compact manifold x 
S^, the S*^ component ?70 of the deformation must be constant. In particular, the 
other two components satisfy 

5{7]lUJ^ + ?72W^) = 

so we conclude that the kernel of L must have dimension bi(AI). 



10 



MICAH WARREN 



4. A politician's optimal transportation problem 

We give the following application of solutions to the above problem. Suppose you 
are the political leader of a compact manifold with non simply connected topology, 
and you are in charge of solving an optimal transportation problem. Due to forces 
in effect before you came into power, you discover that the distributions which need 
to be paired are in fact the same (this fact is known only to your ofhce.) You are 
aware of the unspeakably high political "cost of doing nothing," as well as the clout 
of the Transportation Lobby, so proposing the trivial identity transference plan is 
not an option. Instead, you must come up with a plan which is nontrivial, but does 
not appear arbitrary. Deforming according to a harmonic 1-form as described in 
section 3, you can find a transportation plan which has positive cost, and locally, 
to each supplier in simply connected districts, is given by the cost exponential of a 
potential function, so appears optimal. As long as no one exhibits a set of points 
for which the plan is not cyclically monotone (cf jVi Chapter 5]), it will be believed 
that the plan is a good one . If the deformation is small, such a cycle will involve 
a large set of carefully selected points and will be quite difficult to exhibit. 
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